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We study charmonia correlators at finite temperature. We analyze to what 
extent heavy quarkonia correlators are sensitive to the effect of heavy quark 
transport and whether it is possible to constrain the heavy quark diffusion con- 
stant by lattice calculations. Preliminary lattice calculations of quarkonia cor- 
relators performed on anisotropic lattices show that they are sensitive to the 
effect of heavy quark transport, but much detailed calculations are required to 
constrain the value of the heavy quark diffusion constant. 



1 Introduction 

There are plenty of experimental evidence that strongly interacting matter at high energy density has 
been produced at RHIC [1,2]. One of the most exciting results from RHIC so far is the large azimuthal 
anisotropy of light hadrons with respect to the reaction plane, known as elliptic flow. The observed ellip- 
tic flow is well described by ideal hydrodynamics [3, 4, 5] suggesting early equilibration of the produced 
matter and very short transport mean free path. This interpretation of the experimental data can be chal- 
lenged by measuring elliptic flow of charm and bottom mesons [6, ?, 7]. The first experimental results 
show a non-zero elliptic flow for these heavy mesons. Naively, since the quark mass is significantly 
larger than the temperature of the medium, the mean free pass of heavy mesons is ~ M/T longer than 
the light hadron mean free path. Quantitatively the mean free path is described by the heavy quark dif- 
fusion constant which can be defined through the diffusion equation for the heavy quark number density 
N(x, t), dtN + DV 2 N = 0. If the heavy quark diffusion constant D > 1/T, the predicted heavy quark 
elliptic flow will be too small and in contradiction with current experimental data [8]. 

Kubo formulas relate hydrodynamic transport coefficients to the small frequency behavior of real 
time correlation functions [9, 10]. Correlation functions in real time are in turn related to correlation 
functions in imaginary time by analytic continuation. Karsch and Wyld [11] first attempted to use this 
connection to extract the shear viscosity of QCD from the lattice. More recently, additional attempts 
to extract the shear viscosity [12, 13] and electric conductivity [14] have been made. It turns out that 
Euclidean correlations functions are remarkably insensitive to transport coefficients. For weakly coupled 
field theories this has been discussed by Aarts and Martinez Resco [15]. For this reason, only precise lat- 
tice data and a comprehensive understanding of the different contributions to the Euclidean correlator can 
constrain the transport coefficients. It appeals that heavy quarkonia correlators are the likely candidates 
for meeting this conditions. 

2 Euclidean and real time correlators 

On the lattice we calculate correlation function of local meson operators ( currents ) jfe (x, r) = g(x, r)Thq 
at finite temperature 



with T/j being some combination of the Dirac matrices. This correlation function is related to the real 



time correlation functions D>(x,t,T) = (J h (x,t)J h (0,0)}, D<(x,t,T) = (J h (0,0)J h (x,t)). 
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The most important channels for our further discussion are the pseudo-scalar, = 75 and the 
vector, Th = 7^ channels. In the vector channel the Euclidean correlators are related to density-density 
correlator D^ N = (iV(x, t)N(0, 0)) and current-current correlators Dj} 3 = (J*(x, t) J j (0, 0)), 



G uu (x,r,T) = (J"(x,r)4(0,0)) = —D> N (x., —it, T) , (1) 

G^(x,r,T) = (4(x,r)4(0,0)) = D>} ij (x, -ir,T) . (2) 

Similarly for the pseudo-scalar channel G 5 (x,r,T) = ( j|(x, r) j|(0, 0)) = Df(x, -ir,T) . The 
minus sign in Eq. (1) comes from the relation A = —iA^ between the temporal component of the 
vector in Minkowski space and Euclidean space, in particular x° = —ix% = —it. The spectral function 
is defined through Fourier transform of D£ and Dj^ or equivalently as imaginary mart of the retarded 
correlator x^(k, u) 

D>{\t,u,T) -D<(\t,u,T) 1 

Z7T 7T 

Using the Kubo-Martin Schwinger (KMS) relation £>^(k, t) = D^(k, t+i/T) one discovers the relation 
between the spectral density and the Euclidean correlator, 

/■oo COsh ( OJ I T — T^p ) ) 

Here r is number of zeros in the space-time indexes. 

3 Lattice results on the charmonia correlators and spectral functions 

Charmonia correlators have been studied in lattice QCD and the corresponding spectral functions were 
reconstructed using the Maximal Entropy Method (MEM) [17, 18, 19]. These studies showed that the 
IS states (r] c and J/ip) survive in the plasma up to temperatures as high as 1.6T C . Though it is quite 
difficult to reliably reconstruct the spectral functions, the temperature dependence of the correlators can 
be determined quite precisely [19]. 

We calculated charmonia correlators on quenched anisotropic lattices using the Fermilab formu- 
lation for heavy quarks [16]. Calculation were done at (3 = 6.5 and £ = a s /at = 4, corresponding to 
temporal lattice spacing a^ 1 = 14.12GeV when we set the spatial lattice spacing a s using the Sommer 
scale ro = 0.5fm. We collected about 1000 gauge configurations at each temperature. From Eq. (4) it is 
clear that the temperature dependence of the correlator G(k, r, T) comes from temperature dependence 

of the spectral function and temperature dependence of the kernel K{t,uj,T) = ^ T ~ ■ To sep- 

arate out the trivial temperature dependence due to the kernel K(t, oj, T) following [19] we introduce 
the reconstructed correlator 



G? ec (k,r,T)= / du J p h (k,oj,T = 0)K(T,u J ,T). 
Jo 



(5) 



If the charmonia spectral function do not change across the deconfinement transition temperature T c we 
expect G h /Gf CC ~ 1. In Fig. 1 we show the temperature dependence of G h /G^ cc for pseudo-scalar 
and vector channels at zero spatial momentum k = 0. In the vector channel we show both sum over all 
spatial components J2i G n and the sum over all four components ^2 G^. We see that the temperature 
dependence of the vector and pseudo-scalar correlarors is quite different. For T = 1.5T C we see only 
very small deviations from unity for G h /G^ cc in the pseudo-scalar channel while significant deviations 
are seen in the vector channel. In fact similar temperature dependence of the vector correlator was seen 
in the previous study based on fine isotropic lattices [19, 20]. This is quite unexpected as r\ c and J /if) 
should have similar properties both in the vacuum and in the medium. We will give an explanation for 
this difference in the next section in terms of heavy quark transport. 
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Fig. 1: The ratio G/G rcc for pseudo-scalar (left) and vector (right) channels at k = 0. For the vector case we show both the 
sum over spatial components (filled symbols) and all four components (open symbols). 



4 Spectral functions at low energies and heavy quark transport 

As vector current is a conserved current there should be transport contribution to the corresponding 
spectral function. In general the vector spectral function can be decomposed in terms of transverse 
p T (k, oo, T) and longitudinal jO£,(k, w, T) components. Since the heavy quark mass is much larger 
than the temperature M > Twe can write p ' (k, u,T) = p l( ^ w (k, u,T) + p hi ' gh (k, u, T), where 

p^ h (k,uj,T) contains the resonances and the continuum, and is non-zero for energies lo ~ 2M, and 

L T L T 

Plow (ki 1 ^!^ 1 ) i s tne transport contribution. The simplest way to estimate p^ w (k, lo,T) is to evaluate 
the vector correlator at 1-loop level [21]. In the k — > limit we have p T {<d,uo,T) = p L (0,uj,T) = 
p n {0, u, T), and considering small energies, w«Twe get 

rf ow (0,o;,T) = Xs (T)^uS(u), p? o ° w (0,^,T) = Xl (T) W «( W ). (6) 

Here Xs (T) is the static charm number susceptibility, which in the limit M 3> T is given by Xs (T) = 

12 (^§f^ ^ e~ M / T . Thus at finite temperature we expect that the J2iG n should be enhanced by a 
constant contribution 2>x s {T)T/M relative to its T = value, while the ^^G^ should reduced by 
— X S (T)(1 — 3T/M) (recall Eq. (1)). This is exactly what the lattice data in Fig. 1 show. Furthermore, 
from data on J2i G u and ^ we can estimate that M/T ~ 6 at 1.5T C . The 1-loop result for the 
vector correlator can be also obtained using collision-less Boltzmann equation describing free streaming 
of heavy quarks with no interaction with the plasma [21]. This 1-loop contribution happens to dominate 
the transport part of the Euclidean correlator [21]. 



5 Effective Langevin equation for heavy quark transport 

We have seen that the leading transport contribution to the Euclidean correlator is just a constant and 
corresponds to free streaming of heavy quarks. To get the transport coefficient we need to include the 
effect of heavy quark interactions with the medium. It is very difficult problem in general. Luckily, 
the case of heavy quarks is special since the time scale for diffusion, M/T 2 , is much longer than any 
other time scale in the problem. In terms of the spectral functions, this separation means that transport 
processes contribute at small energy, uj ~ T 2 /M. For this reason we will assume that the Langevin 
equations provide a good macroscopic description of the dynamics of charm quarks [8], 

% = Ji^ = " > «W(f )) = ^5(t - t>) . (7) 

The drag and fluctuation coefficients are related by the fluctuation dissipation relation rj = ^MT- F° r 
time scales which are much larger than 1/rj the heavy quark number density obeys ordinary diffusion 
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Fig. 2: The low energy contribution to the longitudinal spectral function pj^ w (k, u), T) (left) and the low energy contribution to 
the correlator (right). Here uj = uD{M/T) and k = kD^/D/T. 



equation dfN + DV 2 N = 0. The drag coefficient rj can be related to the diffusion coefficient through 
the Einstein relation D = = 

The effective Langevin theory can be derived from kinetic theory in the weak coupling limit [8] 
and probably is adequate for describing heavy quark diffusion even for strongly interacting plasma. The 
Langevin equations make a definite prediction for the retarded correlator \h at small uj and thus for 
the transport part of the spectral functions [21]. The results of calculation for the longitudinal spectral 
functions is shown in Fig. 2. For the case of zero spatial momentum k = we have 

UJ M TT UJ + Tj UJ 

From Eq. (8) it is clear that to calculate the transport coefficient we have to determine the curvature of 
G M (k = 0,r, T) at t = 1/(2T) due to the low energy part of the spectral function pf*. If p\^ w was 
the only contribution to the spectral function and rj = the correlator would be constant. The question 
is how to determine the small curvature in G M (k = 0, r, T), arising from finite value of rj, from the 
curvature arising from the resonance and continuum contributions. This can be done by introducing a 
small chemical potential for the heavy quark, p <C M. Since the transport contribution is proportional to 
Xs, the small chemical potential will enhance the transport by factor of cosh(/u/T) [21]. The small charm 
chemical potential will not affect the resonance and continuum contributions to the spectral function to 
leading order in the heavy quark density, ~ e -( M -v)/T _ Thus we expect that 

6G u = G u ^ r , M ) - G u (r, T, 0) 

Jo m=o smh(uj/{21 )) 

is largely insensitive to the high frequency behavior of the spectral function. In Fig. 2 we show SG U 
quantity for several values of D. The numerical results on the vector correlator show that it can be 
calculated with 0.5% statistical accuracy. Thus if similar numerical accuracy can be achieved for the 
difference 5G n , the curvature and thus the rj, or equivalently D can be estimated in lattice QCD. 
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